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The verification and quantification of experimentally created entanglement by simple measurements, espe- 
cially between distant particles, is an important basic task in quantum processing. When composite systems 
are subjected to local measurements the measurement data will exhibit correlations, whether these systems are 
classical or quantum. Therefore, the observation of correlations in the classical measurement record does not 
automatically imply the presence of quantum correlations in the system under investigation. In this work we 
explore the question of when correlations, or other measurement data, are sufficient to guarantee the existence 
of a certain amount of quantum correlations in the system and when additional information, such as the degree 
of purity of the system, is needed to do so. Various measurement settings are discussed, both numerically and 
analytically. Exact results and lower bounds on the least entanglement consistent with the observations are 
presented. The approach is suitable both for the bi-partite and the multi-partite setting. 
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I. INTRODUCTION 

The theoretical and experimental exploration of entangle- 
ment and in particular its characterization, verification, ma- 
nipulation and quantification are key concerns of quantum in- 
formation science |Q]]. The resource character of entangle- 
ment is most clearly revealed when dealing with situations in 
which a locality constraint is imposed, i.e. when distributing 
the state in such a way that subsequent quantum operations 
can only act on individual constituents supported by classical 
communication. This does not only impose constraints on the 
manipulation and exploitation of entanglement but also on its 
verification. 

In any experiment we will aim to verify the presence of 
entanglement by taking measurements. These measurements 
may either serve to reconstruct the entire state or may only 
collect partial information that is sufficient to reveal the de- 
sired entanglement properties Jj, H, S IH] . Given that a fun- 
damental goal in quantum information science is the creation 
of entanglement between spatially separate locations one is 
often forced to assume that these verification measurements 
are local as well. Generically in such verification experiments 
we will observe correlations in the measurement record. It 
is then a natural question whether these correlations originate 
from quantum correlations in the underlying state or can be 
explained by a classically correlated separable state. Then, if 
there are quantum correlations, one can ask how much quan- 
tum correlations are guaranteed to be there, given the mea- 
surement data. 

Consider as an example a two-qubit system and the mea- 
surement of correlations between Pauli-operators along the z- 
axis, i.e. the quantity 

C zz (pab) = Tr[p AB (o- z ® <t z )] - Tr[p A o- z ] Tr{p B a z ], (1) 

where pa(pb) are the reduced density operators resulting 
from tracing out party B (A) in the original state pab- If 
C zz (pab) — —1 then the measurement outcomes are per- 



fectly anti-correlated and are thus exhibiting very strong, al- 
beit negative, correlations. Do such correlations imply the 
existence of quantum correlations in the underlying quantum 
state? To decide this we must address the following 

Fundamental Question: What is the entanglement content 
of the least entangled quantum state that is compatible with 
the available measurement data ? 

Mathematically, this question is formulated as a minimiza- 
tion problem in which the amount of entanglement in the un- 
derlying quantum state must be minimized subject to the con- 
straints imposed by the measurement data as well as by the 
positivity and normalization of the state The measure- 
ment data will be the expectation values of some observables 
Aj or some non-linear function Fi (p) of the density matrix. 
Then the minimal amount of entanglement -Emm under the 
given constraints is given by 

E nlin = mm{E{p) : Tr[p^] = a.^F^p) = f t } (2) 
p 

where the minimisation domain is the set of states p and E(p) 
is the entanglement measure of choice 12]. Note that this for- 
mulation applies equally to the bi-partite as to the multi-partite 
setting. Note that the importance of the minimization of en- 
tanglement in quantum state reconstruction in quantum infor- 
mation theory was also pointed out in the context of Jaynes' 
principle |9fl. 

The mathematical minimization problem formulated by eq. 
(O may be addressed by techniques from optimization the- 
ory (see e.g. jlOlO . If the constraint are all linear and the 
entanglement quantifier is convex then methods from convex 
optimization theory may be applied. More complicated con- 
straints that are not linear in the density operator (e.g. purity 
measures) can complicate matters considerably. Generally it 
will not be possible to obtain analytic solutions to the opti- 
mization problem and techniques to obtain lower bounds or 
numerical approaches must be used. The analytical and nu- 
merical exploration of these issues will be the main purpose 
of this work. 
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If the optimal state in eq. (O is separable, i.e. E m - m = 0, 
then in reply to our fundamental question we must conclude 
that the available correlations in the measurement record do 
not imply quantum correlations in the underlying quantum 
state. It might be, but need not be entangled. Indeed, in the 
example given in eq. (HJ, the least entangled state compatible 
with the observation C zz (pab) = — 1 is given by 

p = | (|01><01| + |10><10|) (3) 

which is clearly a separable state. Therefore, the observation 
of classical correlations for the measurement along one set of 
directions alone is not sufficient for the verification of entan- 
glement. This well-known observation in quantum informa- 
tion science is particularly relevant in experimental situations 
where only a very restricted set of measurement settings is 
available. 

One way forward consists in measuring additional observ- 
ables. For example, one may consider the measurement of 

Chh(pab) = Tt[pab(th ® an] - Ti[p A an] Tr[p B a ri } (4) 

for all spatial directions n. Observation of perfect anti- 
correlations in all of these measurement records then uniquely 
identifies the sing let state \ip) = 4|(|01) - |10)) as the only 
state compatible with all such measurements. This state car- 
ries one ebit of entanglement. 

In other experimental situations it may be possible to as- 
sert that the state possesses a certain minimal degree of purity 
i33ll . e.g. when decoherence rates, or at least upper bounds for 
it, are known. Let us for example assume that we know not 
only that C zz (pab) = —1 but also that Tr[p^ B ] = 1, i.e. that 
the underlying quantum state is pure. Then again it is straight- 
forward to conclude that the only states compatible with these 
two assumptions are of the form 

|^ = i=(|01> + e ^|10>), (5) 

that is, quantum states with one ebit of entanglement. 

These simple examples serve to make two points. Firstly, 
the simple observation of correlations in measurements along 
a single fixed orientation is not enough to guarantee entangle- 
ment in the underlying quantum state. Secondly, additional 
information, be it correlation measurements along different 
directions or information about the purity of the states, may be 
sufficient to ensure that the correlations found in the classical 
measurement record indeed prove entanglement in the under- 
lying quantum state. Needless to say, in general the situation 
is quite involved as the measurement data may be more varied 
than those in the above examples. It should also be noted that 
the local measurement of the correlation functions mentioned 
above often implies that we possess more information than 
just these correlations. Indeed, we will often possess local 
statistics as well, which in turn can be taken into account when 
answering our fundamental question concerning the minimal 
entanglement compatible with the measurement data. Gener- 
ally, when we are provided with an entangled state, then any 



additional information will make it less and less likely that the 
measurement data is compatible with a separable state. 

Our fundamental question is of particular relevance in ex- 
perimental settings in which it is difficult to perform measure- 
ments for an arbitrarily large number of measurement settings, 
as is required for doing full state tomography. This may be the 
case for example in solid state physics, where it is not always 
straightforward to perform arbitrary measurements. Another 
reason may simply be the existence of constraints on the mea- 
surement time, dictated for example by the stability time of an 
experiment (e.g. in interferometric setups in optics) or by the 
decoherence time (in solid state or other implementations). 

The present work shares some relations with HI , HI and 
lfl9ll where similar questions are developed but where empha- 
sis is placed on observables that are obtained from the theory 
of entanglement witnesses 111311 . Other approaches are consid- 
ered in |3 [H El, [13, Gil]. While these, as well as the present 
work, consider the analysis of a specific state, a somewhat dif- 
ferent approach is taken in 12011 . Here the dynamics of the 
gate used to produce entanglement is considered while mea- 
surements are restricted to a single measurement basis. The 
approach is to make repeated measurements during the gate's 
time evolution. This contrasts with our approach which only 
requires to make measurements on the final state, irrespective 
of the process that created it. 

In this paper we will address our fundamental question for 
systems consisting only of qubits, as this is by far the most 
relevant system from an experimental viewpoint. I should be 
noted however that the approach remains valid unchanged for 
qudits or even infinite dimensional systems. We begin with 
an illustration of the general approach in which correlations 
and purity are quantified by quantum mutual information and 
global entropy, respectively. While these quantities are not 
directly experimentally measurable, they allow for the funda- 
mental question to be most easily answered. Then we consider 
the question for correlations between measurements of Pauli- 
operators along a single axis, e.g. the z-axis. In the process 
we prove an inequality between correlations and purity that is 
valid exactly if a two-qubit state is separable and use it to pro- 
vide necessary and sufficient conditions for entanglement to 
be inferred. Subsequently, we consider correlations along two 
different measurement axes, e.g. x-x correlations and z-z cor- 
relations. Finally we consider the situation in which we take 
into account the local expectation values that are obtained in 
most experiments to sharpen the verification of entanglement. 
We finish with some conclusions. 



II. MUTUAL INFORMATION, ENTROPY AND 
ENTANGLEMENT 

To illustrate the general approach that we are advocating, 
we begin by considering a situation in which the known sys- 
tem properties are the entropy of the state (determining the 
state's purity) and the quantum mutual information (deter- 
mining the state's correlation), and in which the entanglement 
measure of choice is the relative entropy of entanglement 112 ill . 
The reason for this choice is that there exists a very simple re- 
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lationship between these quantities, and the solution of the 
minimization problem eq. (f2]l is immediate. 
The quantum mutual information is given by 

I(pab) = S(p A ) + S(p B ) - S(pab), (6) 

and the relative entropy of entanglement is 111 l23i l24ll 

Er(p) = mm{S(p\\cr) : a separable}. (7) 

For a 2-qubit state, the physically possible values of the pair 
(Iab, Sab) are located in a triangle spanned by the points 
(0,0), (2,0) and (0,2) (see Figure [T). That is, I AB > 0, 
Sab > and Iab + Sab < 2; equality in the latter inequality 
is obtained when both reductions pa and ps are maximally 
mixed. 




FIG. 1: State space in the Sab-Iab plane, depicting the minimal 
possible entanglement Er for every point. In region I there is no 
guaranteed entanglement, while in region II one has at least Er > 

(Iab — Sab)/2. 

The solution to eq. (O is obtained by applying an inequality 
lower bounding the relative entropy of entanglement 112511 and 
showing that equality can be achieved for every pair of values 
of (Iab, Sab)- The inequality is 

Er(pab) > max{S(p A )~S(pAB),S(p B )-S(pAB)} 1 (8) 
which directly implies 

Er(pab) > max(0, \(I( P ab) - S(p AB ))). (9) 

The bound is zero in region I (Sab > Iab) and non-zero 
in region II (Sab < Iab)- Equality in region I is obtained 
by diagonal states; they cover region I completely, and as any 
diagonal state is separable, they have Er = 0. Equality in 
region II is obtained by so-called maximally correlated states, 
which are of the form 



Example 3 of 112411 shows that these states satisfy 

Er(pab) = S(p A ) - S(pab) = S(p B ) ~ S(pab)- (11) 

For any given value of Iab and Sab m region II we can find 
a state of the form eq. (TlOb realizing these values. By eq. (0 
and eq. (TTTb this state realizes the smallest possible value for 
the Er given Iab and Sab- 

The upshot of the results obtained here is that knowledge 
of the two quantities Sab and Iab allows one to have much 
better bounds on the entanglement Er than with just knowl- 
edge of the correlations Iab alone. Indeed, without know- 
ing the purity Sab, one has to assume the worst case, being 
Sab = 2 — Iab, in which case the lower bound on Er is 
given by 



Er > 



0, Iab < 1 

Iab — 1, Iab > !• 



say, 



If, on the other hand, the state is known to be pure, 
(Sab = 0) then the much sharper bound 

Er > Iab I '2 



can be obtained. 

In the rest of the paper we will apply the approach illus- 
trated here for studying the main question eq. (O in the con- 
text of experimentally accessible quantities. In the next Sec- 
tion the measure of correlation will be based on measurements 
along the z-axis. It will turn out that without knowledge of 
the purity one cannot find any lower bound on entanglement 
other than the trivial bound E > 0. Thus, while in the present 
Section one can get some information about the entanglement 
from the quantum mutual information without knowledge of 
the purity, in the next Section knowledge of the purity is ab- 
solutely essential. 



III. PURITY AND CORRELATIONS 

In a number of experimental settings it is not straightfor- 
ward to carry out measurements along arbitrary directions. To 
obtain a measure of correlation in those settings, one can for 
example consider the quantity 

Czz(pab) = TilfiAB^z <S) cr 2 )] - Tt[pa<j z ] Tr[p B <T z ], (12) 

which only requires measurements along the particles' z-axes. 
However, in the previous Section we already alluded to the 
fact that knowledge of this correlation measure alone is not 
sufficient to prove the presence of quantum entanglement. We 
will establish that fact in the present Section. Moreover, we 
will show that if in addition the purity of the state is known, 
as quantified by 



V( P 



AB 



(13) 



Pab = a|00)(00| + 6|00)(11| + fe|ll)(00| + (1 



a)|ll)(ll|. 
(10) 



and provided this purity is large enough, then and only then 
can one infer entanglement from the z-correlation measure. 



Now the question is: how pure does the underlying quan- 
tum state have to be so that \C ZZ \ > indeed implies quantum 
entanglement? Or, more precisely: 

When are all states consistent with given values ofC zz and 
V non-separable, and what is the least entanglement compat- 
ible with these values? 

It turns out that the rigorous analytical answer is surpris- 
ingly involved, largely due to the non-linearity of the con- 
straints involved in the minimization problem, especially if 
one is also interested in the actual amount of entanglement 
that can be guaranteed from such measurements. 

As measure of entanglement we have used the logarithmic 
negativity Ejy, because this is the measure that is most easily 
calculated [28]. The log-negativity is defined as 

E n (pab) :=logTr|/|, 

where V denotes partial transposition w.r.t. subsystem B. 

In Figure |2] we present our numerical results on the small- 
est amount of entanglement compatible with given values of 
purity (see eq. (TT~3T >> and of correlations in the measurement 
record (see eq. (fT2l). This numerical evaluation suggests the 
following: 

• Region I does not allow for any physical states. 

• There is a well-defined central region S that does not 
allow to infer the presence of entanglement as the val- 
ues for purity and correlations can be reproduced by a 
separable state. 

• Only in regions Ha and lib is entanglement guaranteed. 
The minimal value of Em in those regions is given by 

E N >\og 2 {l + y/2{Q-l)+C zz ), (14) 

in Region Ha, and 

E N > \og 2 (C zz + y/2Q - 1) (15) 

in Region lib. Here, Q = Tr[p 2 } = (3T + l)/4. 

One may either calculate both bounds and take the minimum, 
or infer which region one is in via the limits 

i -¥ s « s 5( i+ ( i -¥) 1 )- (16) 

which hold for Region Ha. 

We stress that we do not have a complete proof of these 
statements. They were derived - in a rather laborious way 
- starting from an Ansatz concerning the form of the states 
achieving the bounds. This Ansatz was in turn obtained from a 
combination of Monte-Carlo calculations and inspired guess- 
work. While a proof does not seem forthcoming, the numeri- 
cal evidence for correctness of the Ansatz, and of the bounds 
derived from them, is very convincing. The interested reader 
is advised to contact the authors for further details. 

Analytical proof of boundaries - What we do have been 
able to prove is the analytical form of the boundaries of the S 
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FIG. 2: This plot show the numerical results on the smallest amount 
of entanglement that is compatible with given values of purity (see 
eq. Q~3j) and correlations in the measurement record (see eq. M2V). 
The logarithmic negativity En is shown as a contour plot in function 
of the parameters V and C zz . Four regions can be distinguished. In 
region I no state exists that is compatible with the specified values of 
V and C zz . A large central region, denoted by S, does not allow to 
conclude the presence of entanglement. Regions Ha and lib are the 
only ones where all compatible states have non-zero entanglement 
En- For these regions, the minimal En is given by eqns. d!4t and 
(T5J, respectively. 



region, the region where be separable states. They are given 

by 

^££ < V < i _ ±££ (17) 
3 " " 3 V 

Here, the first inequality defines the boundary with Region I 
while the second one defines the boundary with Region Ha. 

To proceed, we treat boundary I and II separately. For 
boundary I we can simplify the form of p that needs to be con- 
sidered quite significantly. To this end note that correlations 
C zz (p) are unaffected by the transformation 

p^ p = ^ (p + (i ® a z )p{t <g> a z ) + (<j z <g> l)p{a z <g> I) 
+ (<j z ® <r z )p(a z <g> a z )) 

ie, C zz (p) — C zz (p) but at the same time the transformation 
from p to p reduces purity as they correspond to pinchings 
ll3lll . A state p that is invariant under the above maps is di- 
agonal. As these maps are local we find that if E(p) = 0, 
then E(p) = 0. To determine the boundary I let us now 
fix a value for C zz and determine the smallest purity com- 
patible with it. If we have a p with a given purity then by 
the above transformations we can find a p with the same C zz 
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and no larger purity that is diagonal. Therefore it is suffi- 
cient to restrict attention from the outset to diagonal p, i.e. 

p = Diag(a, b, c, 1 — a— b— c). Then we fmdC zz — I — 2b— 2c 
and the purity is given by 



V 



(1 



(18) 



Without restriction of generality we assume C zz > (the case 
C zz < can be treated analogously) and one finds that the 
purity is minimized for b = c. This leaves us with the mini- 
mization of the expression 



IV. CORRELATIONS ALONG DIFFERENT DIRECTIONS 

Let us now move away from the use of non-linear properties 
of the density operator such as purities or entropies and con- 
sider only linear functionals, i.e. expectation values of quan- 
tum mechanical operators that are directly accessible to ex- 
perimental detection. Consider the case when we are given 
the quantities 



C Z2 

C xx 



Tr[(aW 
Tr[(aW 



^)p] 



2b 2 



(l-a-2b) 2 -]) (19) 



for 



1 - 46. (20) 



Then the minimal purity compatible with the given C is then 
found to be 



V 



c 2 

^zz 

3 



(21) 



yielding the boundary confirming eq. (fTTI l. 

Determining the boundary II is more involved and is based 
on the observation that for all separable states p we have rf32ll 

1 



Tr[p 2 ] + -C zz (p) < 1. 



(22) 



We first note that Tr[p 2 ] + \C zz {p) is convex in p. Indeed, a 
short calculation reveals that this expression is equal to 

\Pi k \ 2 + 1 + 2( ^22 - P22 + P§ 3 ~ P3S)- 

As every term is convex in p, the total expression is. There- 
fore, the inequality only has to be checked for the extremal 
points of the set of separable states, i.e. for pure product states. 
This, however, is very easy: for product states, C zz = 0, and 
for pure states Tr[p 2 ] = 1, whence the inequality is satisfied 
with equality. 

Now we note that the separable states p = a\00) (00 1 + (1 — 
a) 1 11} (11 1 saturate the bound (f22l . Rewriting this bound in 
terms of V(p) we find V < 1 — 2C zz /3. This then completes 
the proof for boundary II. 

A lower bound for En - As mentioned above, we have not 
been able to prove our lower bounds (TT4l and ( fT3T > so far. Nev- 
ertheless, inequality (1221 suggests that 

E lower = log^(Tr[p 2 ] + \c zz {p)) (23) 

might be a lower bound on the entanglement in all regions. 
Here we define the function logj(x) := max(0, log 2 (x)); 
that is, loga (x) = for x < 1. We will prove eq. d23l in 
subsection IVIII where a general recipe for the derivation of 
such bounds is presented. 



(note that these are different quantities than the one used in 
the previous Section). In this case it is quite straightforward 
to determine the minimal entanglement compatible with any 
choice of C xx and C zz . To see this we first realize that C xx 
and C zz are invariant under the transformation 



-Q,x,y,z 



(o-j ® a i )p(a l ® <Tj) . 



(24) 



Thus for given C xx and C zz we may restrict attention to states 
of the form 



i+c z 
4 






4 
b 



V%-6 




b 

i-c zz 
4 







1+C„ 



(25) 



Let us now consider the case C xx > and C zz > 0. Any 
other choice can be reduced to this one by application of 1 ® 
a x or 1 (g) a z onto the state. 

The requirements for positivity of p are 6 < (1 — C zz )/A 
and C xx /2 — b < (1 + C 22 )/4. From the first requirement 
follows that 1 6 1 < (l+C zz )/4. Thus any amount of negativity 
of the partial transpose of p must arise from ,014 = C xx /2~ 6. 
As we are looking for the smallest amount of entanglement 
compatible with the choice C xx , C zz > 0, we must minimize 
C xx /2 — 6, i.e. maximize 6. This is achieved by the choice 
6 = (1 — C zz )/A; one checks that this choice satisfies the 
second requirement C xx /2 — 6 < (1 + C zz )/4. Then we find 
E m in{C xx , C zz ) = \ogt{C xx +C ZZ ). For general C xx , C zz 
we find 

E min (C xx ,C„) = log+ (\C XX \ + \C ZZ \). (26) 

This result may easily be generalized to the case of three 
correlations 

C zz = Ttf^Wf )p], 
C xx = Tr[(aW®4 2 ))p], 
C yy = 8(7^ )p], 

for which we find 



Emin = l0gj((l + \C XX \ + \Cyy\ + \C zz \)/2). (27) 
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V. LOCAL STATISTICS FROM CORRELATION 
MEASUREMENTS IMPROVES ENTANGLEMENT 
ESTIMATION 

If the sub-systems for which we would like to verify en- 
tanglement are distant, then any measurement strategy has to 
be composed of local measurements. In this way we can, of 
course, still obtain averages such as (a x (£> a x ) by measuring 
local observables (such as a x ) and use these averages to de- 
termine correlations (such as C xx ). While the assessment of 
entanglement wil primarily depend on the values of these cor- 
relations, it is important to note that these local measurements 
will in addition yield local averages (such as (<r x )), which by 
themselves are not useful to assess entanglement, but when 
taken together with the correlation values represent additional 
knowledge that we can and should take account of. Note that 
the question of the verification of the presence of entangle- 
ment in the particular setting considered in this section has 
been addressed in 01611 . The full analytical treatment of the 
quantification of the least amount of entanglement compatible 
with the measurement data in this setting is quite complicated 
due to the large number of possibilities that are available. In 
the following we will simply present an example to illuminate 
the impact that additional local information may have on the 
question of assessing least entanglement compatible with the 
measurement data. 

Let us reconsider the case in which we employed C zz = 
Tr[(ai 1} <g> a z 2) )p] and C xx = Tc^pP ® a x 2) )p}. In this 
setting we found that E min (C xx , C zz ) = \ogt{C xx + C zz ) 
(eq. (f26|). Let us now investigate what can be gained by tak- 
ing into account knowledge of z\ :— Trfcr^'p] and z 2 := 

(2) 

Tr \a\ >p\\ that is, we determine the minimal amount of entan- 
glement compatible with the information given in C xx , C zz 

and z±, z 2 . 

We can no longer restrict ourselves to states of the form 
d25l l. because Z\ and z 2 are not invariant under transforma- 
tions d24b . The optimal states can now be assumed to possess 
a <7 Z ® a z symmetry. The diagonal elements of the optimal p 
are fully determined by C zz , z\, z 2 and Tr[p] = 1. Employing 
the a z ® a z symmetry of the system the problem can be re- 
duced to a single-parameter minimisation. The optimal states 
turn out to be of the form 



/ a 





/ 


b 


e 





e 


c 





V/ o 





d 



(28) 



and 

< e < Vbc 
< / < Vad. 

Given z\ and z 2 , there are now restrictions on the values of 

C xx and C zz : 

C zz < 1 - \z 2 - Zi\ 

c xx < \^{i + c zz y -{ Zl + z 2 y 
+\^{i-c zz y-{ Zl - Z2 y. 

The negative eigenvalue of the partial transpose of p is given 

by 



A mi „ = - min (1 + C zz - y/(z 1 + z 2 ) 2 + (4e) 2 , 

1 - C„ - VCzi-z 2 ) 2 + (4/) 2 ). 

The log-negativity En is then 

E N = log 2 (l - 2(min(0, A min )). 

To minimise En, we have to maximise A m i n over all allowed 
values of e, which is the range 

max(0, C xx /2 — Vad) < e < mm(C xx /2, \fbc). 

As an example, in Figure [3] we present the difference be- 
tween the minimal compatible entanglement for given (C zz — 
\ — \z\ — z 2 \ = 0.9, C xx , z\ = 0.3, z 2 — 0.2) and the one when 
only (C zz = 0.9, C xx ) are given. For the given value of C zz , 
either b = or c = 0, so that the only allowed value for e is 
e = 0, giving 



A min = -(i-c zz - y/( Zl - z 2 y + (2c xx yy 

While, of course, the parameter range for which physical 
density operators compatible with those data exist is more lim- 
ited in the former case, it is indeed apparent that the knowl- 
edge of C xx and C zz in combination with z\ , z 2 allows us to 
infer a larger amount of entanglement. 

This example highlights the importance of including all 
available information in the entanglement verification as it 
may substantially alter our conclusions. The exact details of 
the procedure will, of course, depend on the concrete situa- 
tion. 



with 



VI. MULTI-PARTITE ENTANGLEMENT 



a = (l + Z! + z 2 + C zz )/4 

b = (1 - zi + z 2 - C zz )/4 

c = (1 + z x - z 2 - C zz )/A 

d = (l-zi-z 2 + C zz )/4 

6 + / = C xx /2, 



Our considerations are not restricted to bi-partite entangle- 
ment. Again, quite general observables may be considered 
but in line with the bi-partite case we illustrate this setting for 
a simple set of observables. Let us consider the expectation 
values C xxx = (<r x <£> a x <£> <j x ), C lzz = (1 ® a z ® a z ) and 
C zz \ — (u z <g) <j z <g) 1 ) . Given the symmetries that leave these 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



c 

xx 

FIG. 3: For given zi := Ti[a z ^ p] and 22 = Tr[cr z 2 ' p] we plot 
the difference between the minimal entanglement (quantified by the 
logarithmic negativity) compatible with the observation of C zz = 
Tr^ai 1 ' ® ai 2) )p] and C xx = Tr^cr^ ® ai 2) )p] and the minimal 
entanglement without the constraints imposed by Z\ and zi- In this 
example, C zz = 0.9. The lower curve is the minimal entanglement 
without knowledge of zi and Z2 (hence the worst case was assumed); 
the upper curve is the minimal entanglement with z\ = 0.3 and 
z 2 = 0.2. 



expectation values invariant, we may restrict attention to den- 
sity operators of the form IB5I1 
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(29) 



In the tripartite setting it is considerably more difficult than in 
the bi-partite setting to define entanglement measures HI . We 
consider two entanglement measures, the relative entropy of 
entanglement and the robustness of enanglement. 

We begin with the relative entropy of entanglement with 
respect to Tri-PPT states, i.e. states that are PPT with respect 
to any of the three possible bi-partite cuts 

E 3 {a) = M{S(a\\p) : pis Tri-PPT}. (30) 

p 

It is helpful to note that it is always sufficient to restrict the 
minimization over p to those states that possess the same local 
symmetries as a I23I l24ll . Thus only states p of the form eq. 
d29l need to be considered. These states all commute with a. 
Thus we are looking for a two-fold minimization 

E min = min{min{S(<r||p) : pis Tri-PPT} : (31) 

a p 

Tr[aAi\ = a.,}. 



We note that states a of the form eq. (|29l are Tri-PPT if and 
only if A = max{|e|, |/|, \g\, \h\} — min{et, b, c, d} < 0. Due 
to unitary invariance of the relative entropy we can apply local 
unitaries to both p and a; one can therefore restrict to non- 
negative real e, /, g and h. Defining m := min{a, b, c, d} we 
obtain the restrictions e, f,g,h< m. 

The expectation values for such states are given by 

C xxx = 2(e + f + g + h) 

C lzz = 2(a-b~c + d) 

C zzl = 2(a + b-c-d) 

1 = 2{a + b + c + d). 

Note that these expectation values lie in the range [—1,1], 

The minimisation over p reduces to a three-parameter min- 
imisation. Let the matrix elements of a and p (in the form 
(|29i >) be denoted a a , a p , etc. The three parameters are a p , b p 
and c p . The other matrix elements of the optimal p are given 
by 

d p = 1 - (a p + b p + c p ) 
hp = win(m p ,(h a /a a )a p ) 
g p = mm(m p ,(g a /b a )b p ) 
f p = min(m p , (/< J /c <T )c p ) 
e p = min(m p , (e ff / 'd a )d p ), 

where m p :— min{a p , b pi c p , d p }. The expression for the rel- 
ative entropy in this optimal state is 

^(crl I p) = H((a a + h a ,a a - K)\\{a p + h p , a p - h p )) + ... 

(32) 

with three additional terms of obvious form. Here, H is the 
classical (Kullback-Leibler) relative entropy between two (un- 
normalised) two-dimensional probability vectors. 

Because of joint convexity of the relative entropy, and con- 
vexity of the feasible set for p, the remaining minimisation 
(over p and a) is a convex one, which means that there can 
only be one local minimum. It can therefore be efficiently 
calculated numerically using, e.g. conjugate gradient meth- 
ods. We have performed numerical calculations based on this 
method, and plotted the results in Figure |4] for the example of 

Cxxx !■ 

Another possible entanglement quantifier is the random ro- 
bustness R(a) i34ll . The random robustness is defined as the 
minimal amount of the maximally mixed state t/d that needs 
to be mixed with a to make the resulting state Tri-PPT. For- 
mally, 

R(a) = min{p : pi / Tr[l] + (1 - p)a is Tri-PPT}. (33) 

p 

We find that 

R(a) = max (o, jj^^j ■ (34) 

Therefore, the minimal robustness under the constraints 

Cxxx = (<J X ® <Jx ® <r x ), Ci zz = (1 ® a z ® a z ) and 
C zz \ = (a z 18) o z (g) 1) is given by 

Rmin = max [0, /Q Am '" ) , (35) 

\ l/0 + A m in/ 
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FIG. 4: For the given value of C xxx = 1, we plot the minimal 
amount of entanglement, as measured by the relative entropy of en- 
tanglement E3 w.r.t. Tri-PPT states, consistent with the observation 

of Cizz and C zz \. 



where 

a r \Cxxx \ — 1 . \Czzl\ + \G\zz\ n1 

A min = max{ - 1 , 0} . (36) 

VII. A GENERAL STRATEGY FOR LOWER BOUNDS ON 
THE NEGATIVITY 

In this Section, we readdress some of the issues of Section 
Hill It is worth noting that the last result obtained there, eq. 
d27l i. could have been obtained from a general strategy to ob- 
tain lower bounds for the minimization problem eq. (|2). 

This can be achieved by using the fact that |jp |ji — 
max|| M || oo=1 Tr[Mp r ] = max|| M || _i Tr[M r /?], where the 
maximization is over Hermitian M ll3lll . Thus we consider 
the problem 

£min = log 2 min{max{Tr[Af r p] : {{M^ = 1} 

p M 

:Tr[pA i ]=a i ,F i (p)<f i },(37) 

where the outer minimization is over positive semidefinite ma- 
trices p (the trace condition for states is included by putting 
Aq = 1, ao = 1), and the inner maximisation is over all Her- 
mitian matrices M. When Fi is a convex function its level sets 
{p : Fi(p) < fi} are convex sets, and we can use the minimax 
equality (see e.g. lllOlO to interchange inner and outer optimi- 
sations, obtaining 



Let us now consider the case that there are no non-linear con- 
straints Fi, then the inner minimization is a semidefinite pro- 
gram (SDP). We now apply Lagrange duality to this mini- 
mization, i.e. we consider the unconstrained minimization of 



the Lagrangian min p Tr[(Af r - J2 v iAi)p] + J2i u i a i over 
all positive semidefinite p > 0, where the Vi are the La- 
grange multipliers. If A/ r — ^ ViAi has negative eigenval- 
ues, the minimum of the Lagrangian will be — oo (by let- 
ting p become arbitrarily large), and will not contribute to 
the outer maximization over M, Thus we can safely require 
AI r — v %Ai > 0, in which case the minimum is obtained 
for p = and equals J2i v i<H- Inserting this we find 

-Emin > log 2 max { max{^ z^a, : ViA t < M r } 

i 

: HMHoo = 1}. (39) 

Because the inner minimization is an SDP, if the problem is 
strictly feasible, i.e. if all inequality constraints can be satis- 
fied with strict inequalities, then we have strong duality ifloll 
and the above step does not weaken the lower bounds. 

Any choice of M and V{ such that M r > J^. v±Ai and 
1 1 M | |oo = 1 now yields a lower bound on E m i n . Indeed, this 
could have been read off immediately from eq. d38l . However, 
as the optimizati on p roblem eq. d39l shows this may be overly 
restrictive. See ll ll [l2tl for lower bounds on other entangle- 
ment measures. 

Applications - In the case of given (C xx , C zz ) with C xx + 
C zz > 1 as discussed in Section ITVl we have Aq = 1, A\ = 
a x o x and A2 = <j z ® <j z , and a,Q = 1, ai = C xx and 
C12 = C zz . In this case we find as optimal M: 

1 1 \ 



\ 1 / 

(which indeed has operator norm 1) and as optimal vi. vq = 
Q, vx = i/ 2 = 1. One checks that M r > A± + A 2 . From this 
we recover again eq. d26l l. 

For the case of given (C xx , C yy , C zz ) we choose i/q = 
vi = v 3 = 1/2, V2 = -1/2, so that J2 l »iM = (1 + <*x ® 
a x — (j y ® (j y + a z ® <t z )/2. Taking M = (1 + <r x ® a x + 
a y ® a y + a z ® cr z )/2 (which has operator norm 1) yields 
M r = '^2 i ViAi, and we recover the exact value found in eq. 

Proof of eq. d23l ) - A similar approach may suggest itself 
for the case concerning purity and correlations discussed in 
section III and will be used to prove the lower bound eq. d23l . 
The constraints are however non-linear. To proceed, we will 
use a kind of linearization procedure. We begin by rewriting 
the quantities C zz (p) and Tr[p 2 ] in terms of expressions linear 
in the tensor product r := p ® p. Taking into account Tr[p] = 
1 we find 

= Tr[rF] (40) 



Z = a z ®a z ®l®l-l®a z ®(T z ®l 

and F is the flip operator that interchanges parties 1, 2 of the 
first copy with parties 3, 4 of the second. The Z presented 



S min = log 2 m A ax{min{Tr[M r p]:Tr[pA i ]=a l ,F l (p)</ l } C zz (p) = Tt[tZ] and Tr[p 2 ] 

: 1 1 I |oo = l}] (38) where Z is the operator 
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here is the simplest one that represents C zz . However, it is 
beneficial to use the symmetrised form Z' = (Z + FZF)/2. 

Let us now address the minimization of | log 2 | \p T <S> p r \ |i 
given constraints on C zz (p) and Tr[/> 2 ]. This problem is lin- 
ear in a = p eg) p and is therefore an SDR Consequen- 
tially, we can apply the above approach. Indeed, let us 
choose M = Diag(llll 1101 1011 1111). Then, clearly, 
M v - Z'/2 - F > and we obtain eq. d23]l as a lower 
bound on the entanglement. This bound is certainly not tight, 
however. Indeed, we could not have expected much more, as 
the extension of the problem to two copies allowed for much 
greater freedom in the matrix AI and, therefore, led us to un- 
derestimate the true value of E min . 



VIII. VERIFICATION OF OTHER PHYSICAL 
PROPERTIES 

In this work we have pointed out that in an experimental 
verification of entanglement we need to search for the least 
entangled state compatible with the measured data. If the state 
so identified is entangled then the experimental data prove 
the presence of entanglement. This approach is not restricted 
to the verification of entanglement. In fact, it applies to any 
physical property that we cannot or chose not to measure di- 
rectly. 

Consider the property II of a quantum system which is 
quantified by IT(p). If we are obtaining experimental data, for 
example quantum mechanical averages of some observables 
Ai, then we need to answer the 

Fundamental Question: What is the least value of II for 
which there is a state that is compatible with the available 
measurement data? 

This smallest value of II is the value to which we have ver- 
ified the presence of IT. Mathematically this may again be 
formulated as a minimization problem in which the property 
IT in the underlying quantum state must be minimized subject 
to the positivity, Hermiticity and normalization and measure- 
ment data obtained as expectation values of observables Ai or 
some non-linear function Fi (p) of the density matrix. Then 
the minimal amount of entanglement E m i u under the given 
constraints is given by 

n min = min{n(p) : Tr[pAi] = a u Fi{p) = fi} (41) 



where the minimisation domain is the set of states p. 

In this more general framework the minimization of entan- 
glement is merely a special case of a general approach to the 
verification of physical properties in experiments. 



IX. SUMMARY AND CONCLUSIONS 

In this work we have addressed the question of when cor- 
relations or other measurement data that have been observed 
in the classical measurement record of a quantum system im- 
ply the existence of quantum correlations in the underlying 
state. The fundamental question in this area may be formu- 
lated as: What is the entanglement content of the least en- 
tangled quantum state that is compatible with the available 
measurement data ? We have formulated this question mathe- 
matically as an optimization problem, discussed it for various 
examples and provided some techniques for obtaining non- 
trivial lower bounds on the minimal entanglement compatible 
with the measurement data. The approach is equally valid in 
the bi-partite and the multi-partite setting and for sub-systems 
of arbitrary dimensionality. We hope that these investigations 
will be helpful in experimental efforts that aim at the creation 
and subsequent unequivocal verification and quantification of 
the generated entanglement. This should, in particular, apply 
to experimental set-ups where for various reasons only a lim- 
ited number of measurement settings is available. 
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